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ABSTRACT 


Large  amplitude  solutions  are  obtained  for  systems  of  semilinear  reaction- 
diffusion  equations  arising  in  mathematical  ecology  which  describe  the  evolu¬ 
tion  of  two  competing  species.  Their  behavior  is  locally  consistent  with  the 
principle  of  competitive  exclusion.  Such  solutions  are  first  obtained  for  a 
special  class  of  steady  state  equations  in  which  the  two  species  are  assumed 
to  be  exactly  equal  competitors?  large  amplitude  patterns  for  generic  classes 
of  equations  are  then  obtained  by  introducing  various  perturbations  in  the 
relative  competitive  strengths  of  the  two  species.  In  particular,  we  obtain 
(1) ,  travelling  wave  solutions  through  constant  perturbations,  and  (2),  stable 
stationary  solutions  through  spatially  inhomogeneous  perturbations. 
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LARGE  AMPLITUDE  PATTERNS  FOR 
TWO  COMPETING  SPECIES 

Robert  A.  Gardner 

0.  Introduction 

The  topic  discussed  here  is  the  existence  of  large  amplitude  solutions  of 
reaction-diffusion  systems  arising  in  mathematical  ecology  which  describe  the 
evolution  of  two  competing  species.  We  seek  solutions  which  are  locally  con¬ 
sistent  with  the  principle  of  competitive  exclusion.  In  particular,  we  obtain 
under  certain  conditions,  stationary  solutions  which  partition  space  into 
regions  in  which  one  species'  density,  u  ,  is  near  its  carrying  capacity  while 
the  density,  v  ,  of  its  competitor  is  small.  At  the  boundary  of  adjacent 
regions,  there  is  a  sharp  "transition  layer"  connecting  these  two  states.  In 
order  to  obtain  such  solutions,  we  must  impose  a  nongeneric  hypothesis  which 
implies  that  the  two  species  are  exactly  equal  competitors.  Such  solutions, 
are  therefore  structurally  unstable.  However,  by  introducing  various  perturba¬ 
tions  into  the  equations  we  are  able  to  obtain  several  patterns  which  are  most 
likely  stable. 

We  first  introduce  a  constant  perturbation  in  the  relative  competitive 
strengths  of  the  two  species,  so  that  one  species  u  ,  is  now  a  slightly  stronger 
competitor  than  the  other,  v  .  In  this  case,  we  obtain  a  travelling  wave  con¬ 
necting  the  rest  state  (u,v)  =  R  at  x  =  -°°,  in  which  u  excludes  v  ,  to  tine 
rest  state  (u,v)  =  p  at  x  =  +“,  in  which  v  excludes  u  .  In  particular, 

(u,v)  is  a  function  of  the  variable  Q  =  x  +  t?t  ,  and  the  wave  moves  from  left 
to  right  with  increasing  t  .  We  do  not  discuss  the  stability  properties  of  such 
solutions  here;  this  will  be  the  subject  of  a  forthcoming  paper. 

We  next  introduce  spatially  inhomogeneous  perturbations.  If,  for  example,  it 

is  assumed  that  u  (resp.  v)  is  the  stronger  competitor  for  x  -  0  ,  (resp.  x 

we  obtain  a  large  amplitude  stationary  solution  connecting  the  state  R  at 
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DAAG29-80-C-0041 ,  This  material  is  based  upon  work  supported  bv  the  National 
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x  =  to  the  state  P  at  x  =  +°°  .  In  this  case,  we  are  able  to  show  that 
such  a  solution,  (or  cline) ,  is  both  structurally  and  exponentially  asymptot¬ 
ically  stable. 

A  number  of  different  techniques  are  employed  to  obtain  the  results  described 
above.  First,  we  use  the  theory  of  isolating  neighborhoods  and  the  generalized 
Morse  index  to  obtain  certain  bounded,  nontrivial  solutions  of  a  broad  class  of 
equations.  However,  such  topological  arguments  yield  only  crude  information 
about  qualitative  behavior.  In  order  to  obtain  finer  information,  we  must  re¬ 
strict  the  discussion  to  a  particular  (Hamiltonian)  class  of  equations.  To  study 
such  systems ,  we  combine  the  above  techniques  with  the  theory  of  bifurcation 
from  simple  eigenvalues.  Finally,  travelling  wave  solutions  and  dines  are 
obtained  by  linearizing  the  appropriate  equations  about  a  large  amplitude  sta¬ 
tionary  solution,  and  by  applying  the  implicit  function  theorem. 
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1.  Formulation  of  the  problem. 

The  equations  to  be  considered  are  of  the  form 


0  = 

u  +  f(u,v) 

(1) 

« 

XX 

0  = 

vxx  +  9(u'v) 

and 

f  Uv 

=  u  +  f(u,v) 

(2) 

J 

; 

1  Vt 

=  vxx  +  g (u,v) 

for  x  <  IR 1 .  We  shall  assume  that  f  and  g  have  the  form  f(u,v)  =  u  M(u,v) 
and  g (u,v)  =  v  N(u,v).  The  functions  M  and  N  are  the  local  growth  rates  of 
u  and  v  respectively.  The  interaction  is  characterized  by  the  following 


hypotheses ; 

(i)  M(u,v)  <_  0  (resp.  N(u,v)  <_  0)  if  u  (resp.  v)  is  sufficiently  large; 

(ii)  M  (u,v)  <  0  and  N  <u,v)  <  0  ,  for  all  u,v  >  0  . 

v  —  u  —  — 

Condition  (i)  is  a  resource  limitation  condition,  whereas  condition  (ii)  charac¬ 
terizes  the  competitive  nature  of  the  interaction.  Competition  equations  have 
been  studied  by  a  number  of  authors;  see  for  example,  Hirsch  and  Smale,  [8], 
Conway  and  Smoller,  [51,  Gardner  [6],  McGehee  and  Armstrong,  [11],  Getz,  [7], 
and  May,  [13].  The  above  hypotheses  include  a  rather  broad  class  of  equations. 
For  definiteness,  we  shall  consider  only  those  fields  (uM,vN)  with  phase 
diagrams  as  in  Figures  1(a)  and  1(b).  We  shall  refer  to  these  as  case  (a)  and 


case  (b) ,  respectively.  Both  cases  satisfy  (i)  and  (ii);  case  (a)  has  been 
studied  by  Conway  and  Smoller,  [5],  and  Gardner,  [6];  case  (b)  has  been  studied 
by  Getz,  [7].  Explicit  examples  of  such  systems  can  be  found  by  letting 


(3a) 

M  =  -v2  +  6j (u-c) (d-u) , 

N  =  -u2  +  <$2  (v-e)  (f-v) 

(3b) 

M  =  2  +  y  -  uv3-u  , 

Y 

N  =  2  +  p  -  u3v  -  v  , 

P 

3' 


in  cases  (a)  and  (b) ,  respectively.  Here  and  f  are  positive 

constants,  and  c,  e,  y,  p  are  constants  such  that  c  <  d  and  e  <  f  . 

It  will  be  convenient  to  denote  (u,v)  by  U  and  (f(U),  g(U))  by 
Solutions  of  (1)  will  also  satisfy  a  first  order  system  of  the  form 
1  U  =  A  =  (a,b) 


(4> 


A  =  -F (U)  , 


F  (U)  . 


•  v 

or  X  =  $ (X) ,  where  X  =  (U,a) ,  and  $  is  the  right  hand  side  of  (4).  If 
P,  Q,  and  R  are  the  rest  points  of  F  indicated  in  Figure  1,  denote  (P,0,0), 
(Q,0,0) ,  and  (R,0,0)  by  P,  Q,  and  R  ,  respectively.  We  shall  ultimately  obtain 
solutions  of  (4)  which  satisfy  the  auxilliary  conditions 
(5)  lim  X(x)  =  R  , 

X  |-x» 


or 


(6)  lim  X (x)  =  R  ,  lim  X(x)  =  P. 

x  -*-<»  X  -»■+«> 

Such  solutions  are  called  homoclinic  and  heteroclinic  orbits ,  respectively . 

We  shall  also  obtain  travelling  wave  solutions  of  (2);  that  is,  solutions  for 
which  U  depends  on  the  single  variable  E  =  x  +  9t.  The  system  (2)  then 
reduces  to  a  system  of  o.d.e. 's  of  the  form 

(7)  U"  +  F(U)  -  eu'  =  0  ; 

(here  "prime"  is  d/dE).  We  shall  seek  solutions  of  (7)  which  satisfy  the 
condition 


(6)  „ 


lim  (U,U*)  =  R  ,  lim  (U,U')  =  P  . 

£  -+  —00  C  -y  +co 


2.  Existence  of  bounded,  large  amplitude  stationary  solutions. 

In  this  section,  we  shall  obtain  bounded,  nontrivial  solutions  of  (4)  which 
have  one  of  the  rest  points  P  or  R  in  their  u>-limit  sets.  The  technique 
is  to  find  an  appropriate  isolating  neighborhood  r,  for  the  flow  of  (4)  so  that 
the  index  of  the  maximal  invariant  set  in  n  is  well  defined.  This  index  is 
easily  computed,  and  it  can  be  used  to  obtain  the  existence  of  the  solutions 


-4- 


described  above.  A  familiarity  with  the  theory  of  isolated  invariant  sets  will 
be  assumed.  A  concise,  nontechnical  description  can  be  found  in  Conley  and 
Smoller,  [4],  A  thorough  development  can  be  found  in  Conley,  [3],  With  the 
exception  of  the  first  two  lemmas ,  the  material  in  subsequent  sections  is  inde¬ 
pendent  of  the  results  of  this  section. 

LEMMA  2.1.  Suppose  that  (U,A)  is  a  solution  of  (4)  such  that  |u(x)|  = 


(u2(x)  +  v2(x) )  ^ 

<  K 

for  all 

X  £  IR  1  . 

Then  there 

exists  L  >  0 

which  depends 

only  on  K  and 

: 

sup 

|U|<K 

i  F  (U)  | 

such  that 

|a(x)  |  <  L 

for  all  x  e 

JR]  . 

Proof .  This  can  be  obtained  from  the  Agmon-Douglis-Nirenberg  estimates  by 
viewing  U  as  the  solution  of  (1)  on  a  unit  interval  I  »  (xQ,x0+l)  with  in¬ 
homogeneous  Dirichlet  conditions.  In  this  manner,  we  obtain  an  H2(I)  bound 
for  U  independently  of  x0  ,  which  in  turn  yields  a  uniform  bound  for  A  on 
I  by  the  Sobolev  lemma.  □ 

Now  let  na  .  *  £a  *  ar  ,  where  a  *  a  or  b  ,  oT  =  {a  e  1R2  :  | A  |  <  L} 

C  f  Li  C  Li  L 

and  ,  a  =  a,b,  are  the  regions  in  the  U-plane  indicated  in  Figures  2(a) 
and  (2 (b) ,  respectively. 


Figure  2 


f  '•  Ij 

LEMMA  2.2.  Suppose  that  L  and  K  =  y/2  c  are  as  in  Lemma  2.1.  Then  nc  L 
is  an  isolating  neighborhood  for  (4)  in  case  (b) ;  g3  .  is  an  isolating 

C  f  Li 

neighborhood  for  (4)  in  case  (a)  provided  that  the  rest  point  S  of_  F  (see 
Figure  2(a))  is  sufficiently  near  the  origin  if  it  lies  in  the  first  quadrant, 
and  that  in  this  case,  the  tangent  at  S  to  M  =  0  (resp.  N  =  0)  is  close 
enough  to  a  vertical  (resp.  horizontal)  line. 

Proof.  We  have  that  Sn"*  .  =  31°  *  a  u  Ea  *  3o  .  a  =  a.b.  Let  (U.A)  be 
-  c,L  c  L  c  L 

the  solution  of  (4)  with  data  (U0,A0)  e  .  If  | A0 |  =  L  ,  then  by 

C  f  Li 

Lemma  2.1,  U (x)  must  exit  from  Ea  _  at  some  value  of  x  .  Now  suppose  that 

c,L 

Uq  e  .  It  is  easily  verified  that  the  vector  field  -F(U)  points  strictly 

Ct  Ct 

out  of  Z  along  3E  .  This  is  obvious  in  case  (b) .  If  in  case  (a)  S  lies 
c  c 

in  the  first  quadrant,  the  additional  assumptions  in  the  statement  of  the  lemma 
guarantee  that  this  will  be  the  case.  The  proof  is  completed  by  expanding  U (x) 
in  a  Taylor  series  about  x  =  0  .  Since  U" (0)  =  -F(Uo>  points  strictly  out  of 

C i  Q 

Zc  ,  it  follows  that  U(x)  is  either  transverse  to  3Ec  or  externally  tangent 

to  3Ea  at  x  =  0  .  □ 

c 

The  relationship  between  case  (a)  and  case  (b)  is  developed  in  the  following 
lemma. 


LEMMA  2.3.  The  flow  of  (4)  ,  case  (a)  in  n_  .  is  related  by  continuation  to 

C  f  Lj 

the  flow  of  (4) ,  case  (b)  in  n  .  . 

- - — - — — — - - —  c,L 

Proof.  We  must  embed  the  two  flows  in  a  family  of  equations  X  =  $>^  (X)  , 

0  <  /  _<  1,  where  <5>q  =  $  (case  (a)),  4>i  =  $  (case  (b) )  ,  and  such  that  there 

exists  an  isolating  neighborhood  for  the  1th  equation  which  is  also  an 

isolc'ing  neighborhood  for  the  ;;th  equation  for  u  sufficiently  near  1  . 

a  b 

Finally,  we  require  that  hq  =  n  and  that  n .  =  n  . 

C  f  L  1  C  f  L 


and  that  n 


Figure  3 


We  begin  by  deforming  the  zero  sets  in  Figure  (2(a)  to  those  in  Figure  3  (a) 

by  pushing  T2  (resp.  T3)  along  the  diagonal  edge  of  E*  until  it  meets 

(resp.  T4) .  Since  u  =  0  (resp.  v  =  0)  is  part  of  the  zero  set  of  f  (resp.  g) 

the  phase  diagram  in  Figure  3(a)  is  topologically  equivalent  to  that  of 

Figure  2(b),  so  that  we  can  deform  Figure  3(a)  to  Figure  3(c)  as  indicated.  It 

is  easily  seen  that  these  deformations  can  be  performed  in  a  manner  such  that 

na  T  is  an  isolating  neighborhood  for  the  flow  defined  by  each  set  of  equations 
C  f  L 


in  the  above  deformation 


Since  the  zero  sets  of  Figure  3(c)  no  longer  include  the  coordinate  axes, 

ci  b 

it  is  easily  seen  that  n  .  can  be  deformed  to  n  through  a  family  of 

c,L  c,L 

isolating  neighborhoods  2  x  o  .  Finally,  since  the  zero  sets  of  Figure  3(d) 

A  Li 

coincide  exactly  with  those  of  Figure  2(b)  in  £b  ,  we  can  deform  the  vector 

c 

field  in  Figure  3(d)  to  the  field  in  Figure  2(b)  ,  in  a  manner  which  leaves  the 

zero  sets  fixed  throughout  the  deformation.  □ 

If  S  (resp.  S,  )  is  the  maximal  invariant  set  contained  in  na  „  (reso. 
a  b  c  ,L 

b 

n  T),  Lemma  2.3  shows  that  the  index  of  S  ,  denoted  by  h(S  ),  is  the  same 

C  f  Li  cl  EL 

as  the  index  of  S,  .  This  index  is  easily  computed  as  follows, 
b 

LEMMA  2.4.  The  index  of  is  the  pointed  2-sphere  (denoted  by  Z2). 

Proof.  It  is  easily  seen  that  can  be  continued  to  the  single  hyperbolic 

rest  point  p  by  pulling  the  curves  M  =  0  and  N  =  0,  (case  (b) ) ,  apart  in 

a  manner  such  that  the  new  vector  field  equals  F  near  P  and  3Zb  and 

admits  the  unique  rest  point  P  in  £b  .  It  is  obvious  that  {£>}  is  the  onlv 

c 

k 

complete  solution  of  the  new  equations  in  n  „  .  Hence  the  index  of  S,  equals 

C  ,L  b 

the  index  of  {p},  considered  as  an  isolated  invariant  set. 

The  rest  point  P  of  F  is  a  node,  so  that  the  eigenvalues  X  and 
of  dF  at  P  are  both  negative.  A  simple  computation  shows  that  the  eigenvalue 

of  dO  at  p  are  ±/-V  ,  i  =  1,2.  it  follows  that  the  index  of  {p}  is  , 

since  the  linearized  equations  are  nondegenerate  and  have  two  positive  eigenvalue 

We  now  apply  this  result  to  obtain  some  information  about  the  sets  S  and 

a 

sb  .  Let 

Me<X0)  =  {X  e  1R4  ;  0  <  |x  -  X0|  <  el. 

THEOREM  2.5.  For  every  e  >  0  ,  we  have  that  either  S  r.  M  (p)  f  0  or  that 
S  n  M  (R)  =  (6,  where  S  =  S  ,S.  . 

E  a 


-8- 


Proof.  If  this  was  not  the  case,  Sj  =  S  \  {P,R}  would  be  an  isolated  in¬ 
variant  set;  by  an  easy  computation,  (see  Conley,  [ 3] ) ,  we  have  that 
(8)  h (S)  s  h (?)  v  h(R>  v  h (s3 )  , 

where  "v"  is  the  topological  sum  of  pointed  spaces  and  —  denotes  a  homotopy 
equivalence  of  such  spaces.  We  obtain  a  contradiction  by  computing  the  second 
homotopy  groups  of  the  left  and  right  hand  sides  of  (8).  Since  homotopy  equiv¬ 
alence  induce  isomorphisms  on  the  hcmotopy  groups,  by  Lemma  2.4  and  Maunder,  [12, 
Theorem  7.4.6],  we  have  that 

Z  as  tt2  (E2)  —  tt2  (h(P)  v  h(R)  v  htS^). 

However,  for  two  pointed  spaces  X  and  Y  we  have  that  by  Maunder,  [12,  Theorem 


7.2.20], 


tt2(X  v  y)  a  tt2  (X)  ©  rr2  (Y)  ©  tt3  (X  x  y ,  x  v  y)  , 


where  n  (Z  Z  )  is  the  third  relative  homotopy  group  of  the  pair  (Z  ,Z  ).  Hence 
3  l  f  ^  12 

we  have  that 

Z  a  7r„(h(p)  V  h(R))  0  tt  (h(S  ))  ®  G  , 

2  2  1 

where  G  is  tt3  of  the  pair 

(h (p)  v  h(R)  x  h(Si),  h (P)  v  h(R)  v  h (S j ) ) . 

However,  by  the  proof  of  Lemma  2.4  we  have  that  E2  s  h(P)  s  h(R),  and  therefore. 


tt 2  (h (P)  V  h(R))  s  tt2(E2  V  E2)  as  z  ®  z  , 

by  Maunder,  [12,  p.  297],  Thus  2  ^2  ©2  ©  tt  2  (h  ( ) )  ®  G,  yielding  a  contra¬ 
diction,  r 


COROLLARY  2.6.  There  exists  a  nontrivial  solution  of  (4)  in  S  ,  a  =  a.b. 

■  1  ■ '  ■  '  ■  1  ■  11  -  *  ■  ■  a 

whose  us-limit  set  consists  of  one  of  the  rest  points  P  or  R  . 

Proof .  This  is  an  immediate  consequence  of  Theorem  2.5  and  the  Hartman  lineari¬ 


zation  theorem.  For  brevity's  sake,  we  omit  the  details 


3.  Hamiltonian  systems . 

We  shall  now  assume  that  there  exists  a  function  H(U)  such  that  VH  =  F  , 
so  that  if  E(U,A)  =  |a]  '  +  2h(U),  then  E  is  constant  on  solutions  of  (4). 

This  will  be  the  case  if  M  and  N  are  chosen  as  in  (3a)  or  (3b).  (More 
generally,  if  the  zero  sets  of  M  and  N  are  v  =  k(u)  and  u  =  £(v), 
respectively,  let  F  =  (u(-v^  +  k(u)2),  v(-u2  +  2(v)2)).  Then  F  is  a  gradi¬ 
ent  field  whose  components  have  the  same  zero  sets  as  those  of  F) .  The  solu¬ 
tions  ,  now  constrained  to  lie  in  a  three  dimensional  manifold,  have  less  room 
to  move  around,  and  substantially  more  can  be  proved  about  their  behavior.  We 
shall  assume  that  F  is  as  in  case  (b) ,  though  the  following  results  also 
hold  for  certain  systems  of  type  (a) ;  this  is  discussed  at  the  end  of  this 
section. 

Hamiltonian  systems  in  this  context  do  not  occur  generically,  and  their  use 
in  mathematical  ecology,  where  the  equations  are  only  crudely  known,  is  rather 
artificial.  However,  the  assumption  of  such  additional  structure  is  reasonable, 
since  the  intention  here  is  to  prove  complexity,  rather  than  simplicity  of 
behavior.  More  precisely,  the  sign  of  the  quantity  H(P)  -  H (R)  determines 
which  of  the  two  species  is  the  stronger  competitor;  we  will  show  later  that 
certain  qualitative  properties  of  solutions  when  this  quantity  is  zero  persist 
under  various  perturbations  of  the  nonlinear  term.  Furthermore,  in  a  forthcoming 
paper,  we  shall  obtain  travelling  wave  solutions  of  quite  general  classes  of 
equations  by  constructing  a  homotopy  back  to  a  Hamiltonian  system  for  which 
H  (P)  =  H  ( R)  . 

The  strategy  is  to  obtain  small  amplitude  periodic  solutions  (from  bifurca¬ 
tion  theory)  which  have  values  in  .  These  branches  lie  in  globally  de- 

c ,  L 

fined  continua  of  solutions  which  also  have  values  in  .  This  fact,  com- 

c  ,L 

bined  with  the  Hamiltonian  structure,  can  be  used  to  obtain  a  precise  description 
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of  large  amplitude  periodic  solutions,  which  in  turn  can  be  used  to  obtain 
heteroclinic  or  homoclinic  orbits  depending  on  H(P)  -  H(R). 

To  obtain  small  amplitude  periodic  solutions,  we  consider  the  equations 

(9)  X“2U  +  F(U)  =  0  ,  U  (-1)  =  U  (1)  =  0  , 

XX  XX 

where  X  is  a  large  positive  parameter.  If  U(x)  is  a  solution  of  (9),  then 
U(± A  ^x)  are  solutions  of (1)  on  the  interval  -X  <  x  <  X  ,  with  homogeneous 
Neumann  conditions  at  the  boundary.  By  translating  and  "piecing  together"  such 
solutions  in  an  obvious  manner,  we  obtain  a  2X-periodic  solution  of  (1),  (or 
of  (4)),  on  1R1  .  (This  procedure  yields  a  genuine  solution  of  (4),  as  can 
be  seen  by  applying  the  uniqueness  theorem  for  o.d.e.'s  at  points  where  two 
solutions  are  pasted  together) . 

Nontrivial  solutions  of  (9)  near  Q  ,  (the  saddle  in  Figure  1  (b) ) ,  are 
obtained  as  follows.  The  Jacobean  of  F  at  Q  has  one  positive  and  one  nega¬ 
tive  eigenvalue,  p  and  n  respectively.  The  bifurcation  points  of  (9)  will 
occur  at  values  of  X  =  X^  for  which  -n2ir2/4A2  +  p  *  0  .  If  ej  is  a  unit 
is  a  unit  vector  in  IR2  in  the  direction  of  the  unstable  manifold  at  Q  , 

(for  the  equations  U  =  F(U)),  we  obtain  solutions  of  (9)  with  A  *  Xr  +  0( y2) 
of  the  form 

(10)  U  (x,  y)  =  Q  +  yej  sin(nirx/2)  +  0(y2), 

n 

where  y  is  a  small  parameter.  Let  V^(x,y)  be  the  2X (y)-periodic  solution 

of  (1)  obtained  from  Un(x,y)  as  above. 

Now  let  X  be  the  space  of  pairs  of  bounded  continuous  functions  on  (-1,1), 

and  let  £  =  X  x  b+  .  Then  by  Rabinowitz's  theorem,  [18,  Theorem  1.3], 

(U  (x,y),  X  (y) )  lies  in  a  globally  defined  continuum  of  solutions  U ^  £ £,  in 

the  sense  that  either  "loops  back"  to  the  trivial  branch  Q  *  at  some 

value  of  X  =  X  where  m  /  n  ,  or  U  meets  »  in  £  .  We  also  obtain  from 

in  n 
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U  in  the  manner  described  above  a  continuum  V  of  periodic  solutions  on 
n  n 

1R1  containing  (V  (x,y) , A  (y) ) .  It  is  easily  seen  that  if  (V,A)  e  1/  ,  then 

(V(x) ,V' (x))  e  for  all  x  .  From  the  expansion  (10),  this  clearly  holds 

C ,  L 

for  (V,A)  near  (Q,A  ).  If  the  assertion  fails  for  some  (V,A)  e  V  ,  we  can 

n  n 

connect  (V.A)  to  (Q,A  )  with  a  continuous  curve  in  1/  .  We  can  then  use 

n  n 

the  continuity  of  the  V-component  along  this  curve,  our  choice  of  c  and  L  , 
and  Lemmas  2.1  and  2.2,  to  obtain  a  contradiction. 

Rabinowitz  also  characterizes  the  nodal  properties  of  solutions  of  nonlinear 
Sturm-Liouville  problems  which  lie  on  continua  which  bifurcate  from  the  trivial 
solution.  However,  his  argument  doesn't  have  an  obvious  generalization  to 
systems.  We  shall  prove  an  analog  of  this  result  for  systems  in  the  special 
case  under  consideration.  To  do  so,  we  shall  need  the  following  lemma. 


LEMMA  3.1 


Suppose  that  (U,A)  is  a  non-constant  solution  of  (4) ,  that  for 


some  x  =  x0,  U(xQ)  does  not  lie  in  the  shaded  region  in  Figure  4(a),  and  that 
A(x  )  lies  in  either  the  first  or  the  third  quadrants  of  the  A-plane,  (including 
a  =  0  and  b  =*  0).  Then  (U/a)  is  not  a  bounded  solution  of  (4)  with  values 


Proof .  Suppose  that  U(xQ)  is  as  above  at  x  =  xQ.  Without  loss  of  generality, 

we  can  assume  that  U(xQ)  lies  above  the  shaded  region;  (it  is  easily  seen  that 

if  U(x0)  lies  in  the  extreme  upper  left  -  or  lower  right  -  hand  corners  of  i*3  f 

u  c 

the  lemma  is  true  regardless  of  a(xq)).  We  can  also  assume  that  A(x0)  points 
into  the  first  quadrant;  (replace  U(x)  with  U(-x)).  A  moment's  reflection 
shows  that  we  can  actually  assume  that  a(xq)  points  strictly  into  the  first 
quadrant. 

It  clearly  suffices  to  show  that  U(x)  cannot  cross  the  vertical  or  hori¬ 
zontal  lines  through  U(x0)  for  any  x  >  xfi  .  Suppose,  for  example,  that  U(xJ 
lies  on  the  horizontal  line  for  some  x^  >  x^  ,  By  the  mean  value  theorem,  there 
must  exist  x2  between  xQ  and  x1  such  that  a(x2)  points  strictly  into 
the  second  quadrant  of  the  A— plane .  However,  the  mean  value  theorem  applied  to 
A(x0)  and  A(x2)  yields  a  value  x3  between  xQ  and  x2  such  that  -F(U(x.)) 
points  strictly  into  the  half  plane  {b  <  0}  .  However,  for  U  above  the 
shaded  region,  -F(U)  points  into  the  first  quadrant,  yielding  a  contradiction.  — 
LEMMA  3.2.  Let  L  be  the  straight  line  in  the  U-plane  passing  through  Q 
and  through  the  origin.  Then  any  non-constant  element  of  l!j 
crosses  L  exactly  once  for  x  e  [-1,1], 

Proof .  From  the  expansion  (10) ,  the  assertion  is  clearly  true  for  solutions 
which  are  close  to  (Q,Aj).  Suppose  that  the  lemma  is  false  for  some  (U,\)  U.  . 
We  connect  (U,A)  to  with  a  continuous  curve  in  Uj.  By  standard 
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regularity  theorems,  we  see  that  U*  varies  continuously  (in  the  X-topology) 

along  this  curve.  Hence  there  must  exist  an  element  (U,A)  along  this  branch 

which  intersects  L  tangentially  at  least  once.  If  V  is  the  2A-periodic 

solution  corresponding  to  U  ,  then  a  similar  statement  holds  for  V  at  some 

value  of  x  =  x  .  However,  (V(x),V*(x))  satisfies  the  hypotheses  of  Lemma  3.1, 

so  that  (V,V')  cannot  be  a  bounded  solution  of  (4)  with  values  in  n  , 

c  ,L 

yielding  a  contradiction.  □ 

LEMMA  3.3.  The  continuum  U cannot  loop  back  to  the  trivial  branch  at  any 
finite  bifurcation  point  (Q,Am),  A  A  j  >  so  that  llj  must  meet  ®  in  £  . 
Moreover,  the  projection  of  Uj  on  IR+  contains  an  interval  of  the  form 

t^,”)  • 

Proof.  The  first  assertion  follows  from  Lemma  3.2  and  the  explicit  expression 

(10)  for  solutions  near  (Q,A  ),  m  >  1  .  The  second  assertion  follows  from 

m 

Rabinowitz's  theorem.  The  last  assertion  follows  from  the  fact  that  we  have 
a  bound  on  the  projection  of  U on  X  ,  namely,  |u(x) |  <  /l  c  for 
(U,A)  €  Uj  .  □ 

LEMMA  3.4.  Suppose  that  (U,A)  e  U ^  and  that  U  is  not  a  constant.  Then 
(i)  U(il)  lies  in  the  interior  of  the  shaded  region  in  Figure  4(a),  and 

(11)  the  components  of  U  are  strictly  monotone  on  (-1,1). 

Proof.  The  first  assertion  follows  from  the  fact  that  u  (-1)  =  U  (1)  =  0  . 
-  x  x  ' 

and  from  Lemma  3.1. 

We  now  prove  (ii).  For  definiteness,  suppose  that  u(l)  lies  in  the  upper 
shaded  region  of  Figure  4(a),  so  that  by  Lemma  3.2,  U(-l)  must  lie  in  the 

lower  shaded  region.  If  U  =  (u,v) ,  let  xQ  be  the  value  of  x  closest  to 
but  different  from  -1  such  that  u' (x)  =0  or  v' (x)  =  0.  Suppose  that 
u' (xQ)  =  0  and  that  U(Xq)  lies  on  or  below  L  .  (If  U(x0)  lies  above  L, 
then  let  x;  be  the  value  of  x  closest  to  but  different  from  +1  such  that 
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u  or  v  has  a  critical  point;  then  U(xj)  lies  above  L  ,  and  the  following 
argument  can  be  applied  to  U  at  x^). 

Since  U(-l)  lies  in  the  lower  shaded  region  of  Figure  4(a),  we  have  that 
u  (-1)  <  0  and  that  v  (-1)  >  0  ,  so  that  u  has  a  local  maximum  and  v 

XX  XX 

has  a  local  minimum  at  x  =  -1.  Since  u  has  no  critical  points  between  -1 
and  x  ,  xQ  must  either  be  a  point  of  inflection  or  a  local  minimum  of  u  ; 
that  is,  uxx(xQ)  =  -f (U(xQ))  0.  However,  the  region  on  or  below  L  in  which 

-f (U)  >_  0  consists  of  points  that  lie  above  the  lower  shaded  region  in  Figure 
4(a).  Moreover,  U' (xQ)  =  (0,v'(Xg))  is  either  vertical  or  trivial.  If  V  is 
the  periodic  solution  on  IR  1  corresponding  to  U  ,  we  have  for  some  x  =  x 
that  (V(Jf)  ,V' (5c) )  satisfies  the  hypotheses  of  Lemma  3.1.  However,  (V,V')  can¬ 
not  exit  from  r|b  ;  we  have  therefore  obtained  a  contradiction.  □ 

c,L 

THEOREM  3.5.  Suppose  that  E(P)  =  E (R) .  Then  there  exists  a  heteroclinic  orbit 
of  (4)  which  satisfies  (6). 

Proof.  Select  a  sequence  (U  ,1  )  e  U,  such  that  lim  l  =  “.  Let  x  e  (-1,1) 

n  n  1  n  n 

n-x» 

be  the  unique  value  of  x  for  which  Un  (x)  e  L  .  By  passing  to  a  subsequence 

if  necessary,  we  can  assume  that  lim  x  =  x„  e  [-1,1],  that  x  <0  for 

no  n  — 

n  =  0,1,2,...  ,  and  that  U  (1)  lies  above  L  for  all  n  .  If  either  of  the 

n 

last  two  statements  is  reversed,  the  proof  is  similar. 

We  first  show  that  Un(l)  cannot  have  a  subsequence  which  approaches  Q  . 

Assume  that  this  is  not  the  case  and  let  the  subsequence  again  be  denoted  by 

U^d).  Let  vn(x)  be  the  2 J^-periodic  solution  corresponding  to  Un  .  Then 

(11)  E (V  (*.  )  ,0)  =  E (V  (-£  )  ,0) 

n  n  n  n 

for  all  n  ,  It  is  easily  seen  that  for  any  W  ^  Q  in  the  shaded  region  of 

Figure  4(a),  H(w)  >  H (Q) .  Thus  vn^~s-n^  must  also  approach  Q  ,  since  if  any 

subsequence  approached  a  point  W  as  above,  (11)  would  be  violated  for  suffi¬ 
ciently  large  n  .  Since  the  components  of  Un  are  monotone  on  (-1,1),  the 

components  of  V  are  monotone  on  (-2.  .1  ).  Thus  V  converges  uniformly  to 
n  n  n  n 

Q  on  IP 1  . 


We  obtain  a  contradiction  as  follows.  Translate  V  (x)  to  2  (x)  = 

n  n 

Vn(x  +  *„XJ  so  that  2  (0)  e  L  ,  and  if  p  =  «.  (1-x)  >«.  ,  then  Z  (x) 

11  n  ji  n  n  n  n  —  n  n 

lies  above  L  for  all  x  ,  0  <  x  <  p  .  Let  Z  »  Z  -  0  ,  and  let  d  = 

dist  (Z  (p  ),Q).  By  the  monotonicity  of  the  components  of  Z  on  tO.p  )  we 
*  11  n  n 

have  that  Z^(x)  lies  within  a  distance  of  d^  from  L  ,  for  x  in  this  inter 

val.  Now  if  M  =  dF  at  U  =  Q  ,  we  have  that  for  0  <  x  <  p  , 

—  —  n 

(12)  z"  +  MZ  =  0{l2)  =  0(d2). 

n  n  n  n 

Now  let  et  be  the  unit  vector  in  IR2  tangent  to  the  unstable  manifold  of  f 

at  Q  ;  we  have  that  e^  points  strictly  into  the  second  quadrant  and  that 
t  t 

elM  *  P®i  »  where  p  is  the  positive  eigenvalue  of  M  .  Now  let  E(x)  = 
el  sin  (/px) ;  multiply  (12)  by  E  and  integrate  from  x  =  0  to  x  =  ir//p  to 
obtain 

(13)  -e*  •  (Zn(7r/^p)  -  Zn(0))  =  0(d2)  . 


Note  that  since  x  <  0,  we  have  that  <  P  for  large  n  .  Let  ~  be 

n  n  n 

the  set  of  points  lying  above  L  and  at  a  distance  of  not  more  than  ->d  from 

n 

L,  where  a  e  i0,$)  is  a  constant  to  be  chosen  as  follows.  Since  the  tangents 


at  Q  to  M  =  0  and  N  =  0  both  have  negative  slope,  the  line  connecting 
a  point  2^  to  Gn  can  be  made  to  have  negative  slope  by  choosing  cx  suffi¬ 
ciently  small,  (in  a  manner  depending  only  on  the  slopes  of  the  above  two  tan¬ 
gents),  where  Z^  is  any  point  in  the  upper  shaded  region  in  Figure  4(a)  at 
a  distance  of  d^  from  L  ,  and  is  any  point  in  the  upper  shaded  region 

lying  in  r  .  Thus  if  0  is  the  angle  between  2  -  G  and  e,,  then 

cosB^  is  bounded  away  from  zero  independently  of  n  . 

Suppose  that  a  subsequence  of  2  (ir/^p)  lies  outside  r  .  Let  the  subse- 

n  n 

quence  again  be  denoted  by  Z^fir/i/p).  By  the  monotonicity  assertion  of  Lemma 

3.4  the  vector  Z^-n/v'p)  -  Z^(0)  points  into  the  second  quadrant  whereas  e, 

points  strictly  into  the  second  quadrant.  Thus  if  0^  is  the  angle  between 

these  two  vectors,  cos@  is  bounded  away  from  zero  for  all  n  .  Hence  the 

n 

left  hand  side  of  (13)  is  0 (d  ) ,  yielding  a  contradiction. 

n 

We  must  therefore  have  that  2  (ir/i/p)  lies  in  r  for  all  n  ,  and  by 

monotonicity,  a  similar  statement  holds  for  Z^(x)  for  ell  x  €  [0,v/Sp\.  If 

2  (x)  lies  inside  the  upper  shaded  region,  connect  Z  (x)  to  2  =  2  (p  ) 

n  n  n  n  n 

with  a  straight  line,  and  let  G  be  the  point  of  intersection  of  this  line 

n 

with  the  left  hand  boundary  of  T  .  If  2  (x)  doesn't  lie  in  the  upper  shaded 

n  n 

region,  connect  2  (x)  to  a  point  A  on  the  closer  of  the  two  curves  M  =  0 
n  n 

or  N  =  0  with  a  line  parallel  to  L  ,  (and  hence,  in  r  ).  Now  connect  a 

n  n 

to  2  *  2 (p  )  with  a  straight  line  and  let  G  be  as  above.  Since  H  =  f 

n  n  n  3  n  u 

end  Hv  *  g  are  both  negative  (resp.  positive)  above  (resp.  below)  the  upper 

shaded  region,  and  since  <  0  and  >  0  inside  this  region,  we  see  from 

our  choice  of  a,  that  in  both  of  the  above  two  cases,  H(Z  (x) )  <  H(G  )  for 

n  n 

all  x  e  [0,it/</p].  Thus 

i  ,2 

Z* (x) |  -  2H(Z  )  -  2H(Z  (x)) 

n  n 

>  2H (Z  )  -  H(G  ) 
n  n 

-  2VH (U*)  •  (Z  -  G  )  , 
n  n  n 
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where  U*  is  a  point  on  the  line  connecting  to  Z^  .  Since  Q  is  a 

nondegenerate  rest  point  of  dF  ,  we  have  that  there  exists  a  constant  K  >  0 
such  that 

| VH (U)  j  =  |f(U)  |  >_  K|u  -  Qj 

for  U  near  Q  .  Moreover,  F(U*)  points  into  the  second  quadrant,  and  by 

n 

our  choice  of  a  ,  Z^  -  G^  points  strictly  into  the  second  quadrant  in  a  manner 

such  that  if  0  is  the  angle  between  F(U*)  and  Z  -  G  then  cosQ  is  bounded 
n  n  n  n  n 

away  from  zero  for  all  n  .  Hence  there  exists  a  constant  c  >  0  such  that  for 

0  <_  x  <_  Ti/t/p  we  have  that 

|z* (x) I2  >  cd2  , 

1  n  1  —  n 

since  |u*  -  o|  >  ad  and  Iz  -  G  |  >  (1  -  a)d  .  Moreover,  by  the  monotonicity 
1  n  * 1  —  n  1  n  n1—  n 

assertion  of  Lemma  3.4,  we  have  that  Z* (x)  points  into  the  second  quadrant,  so 

n 

that  if  0  is  the  angle  between  Z* (x)  and  e,  ,  we  have  that  cos0  is 

n  n  1  n 

bounded  away  from  zero  for  all  n  and  x  e  [0,7r/v/p].  Hence,  the  left  hand  side 

of  (13)  is  once  again  0(dn),  yielding  a  contradiction.  We  must  therefore  have 

that  |v  (S,  )  -  o|  >  e  for  all  n  and  for  some  e  >  0  . 

We  now  obtain  the  existence  of  a  heteroclinic  orbit  of  (4)  connecting  P 

and  R  .  By  passing  to  a  subsequence  if  necessary,  we  can  assume  that  Z^ (0) 

converges  to  a  limit  V  ,  and  that  Z^(0)  converges  to  a  limit  A  .  Let 

(V(x),A(x))  be  the  solution  of  (4)  with  data  (V,A) ,  and  let  J  =  [0,x],  where 

x  >  0  is  arbitrary.  Since  x  <  0  we  have  that  the  components  of  Z  (x)  must 

n  —  n 

be  monotone  on  J  for  large  n  ,  and  since  Z^  uniformly  approximates  V(x) 
on  J  ,  the  components  of  V(x)  must  be  monotone  for  all  x  >  0  .  Hence 
lim  v(x)  =  V  exists  and  is  finite.  Thus  the  solution  (Vm(x) jA^fx) )  through 

X  w 

any  point  in  the  w-limit  set  of  (V(x),A(x))  must  have  constant  V-component; 
thus  A  (x)  :  0  and  F(V  )  =  0.  Hence  lim  A(x)  =0  and  V  is  a  rest  point 

oo  00  00 

x-*» 

of  F  .  It  follows  from  the  preceding  result  that  ^  Q.  The  only  other 
possibility  is  that  =  P. 
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We  finish  the  proof  by  showing  that  lim  (V(x),A(x))  =  R  .  Let  o  =  (-1-x  ) 

n  n 

x  ► 

and  suppose  that  some  subsequence  of  Z  (o  )  converges  to  a  limit  R,  4  R  . 

n  n  l 

From  Lemma  3.4,  we  see  that  R  must  lie  in  the  lower  shaded  region  of  Figure 

4(a).  since  H  (R  )  <  H (R)  and  E  (Z  (o  )  ,0)  =  E(Z  (p  )  ,0)  -*■  E(P)  =  E(r),  we 
i  n  n  n  n 

obtain  a  contradiction  for  sufficiently  large  n  .  Thus  we  have  that 

lim  Z  (o  )  =  R.  Now  suppose  that  lim  o  >  and  select  a  subsequence  (again 
n  n  n  ~  n 

denoted  by  a  )  such  that  lim  a  =  a„  >  -°°.  Since  (Z  (x) ,  Z'(x))  converges 
n  no  n  n 

n 

uniformly  to  (V(x),A(x))  on  [a0,O],  we  must  have  that  V(aQ)  =  R  and  that 

A(o0>  =  0  ,  so  that  (V (x)  ,A  (x) )  =  R  ,  contradiction.  Thus  lim  a  =  -»  ,  and 

the  argument  of  the  preceding  paragraph  can  now  be  applied  to  obtain 

lim  (V  (x)  ,A  (x) )  =  R  .  □ 

x->— «> 

THEOREM  3.6.  Suppose  that  E(p)  >  E(r)  .  Then  there  exists  a  homoclinic  orbit 
of  (4)  satisfying  (5)  . 

Since  the  details  of  the  proof  are  similar  to  those  of  Theorem  3.5,  they  will 
be  omitted. 

The  above  arguments  can  be  applied  verbatim  to  systems  of  type  (a)  provided 
that  the  zero  sets  of  M  and  N  are  monotone  decreasing  with  increasing  u  . 
This  will  be  the  case  if  the  smallest  roots  of  the  two  parabolas  are  sufficiently 
negative,  see  Figure  6.  In  a  forthcoming  paper,  we  shall  obtain  travelling  wave 
solutions  of  systems  for  which  (uM,vN)  has  such  qualitative  properties  and  for 
which  (uM, vN)  is  not  necessarily  a  gradient. 


v 


We  finally  remark  that  the  method  of  proof  of  Theorem  3.5  shows  that  every 
branch  on  the  continuum  \i j  "approaches"  a  homoclinic  or  heteroclinic  orbit, 
in  the  sense  described  above.  Thus  the  steady  state  equations  (1)  admit  the 
familiar  "teardrop”  of  bounded  solutions  which  arises  in  the  case  of  a  scalar 
equation. 

4.  Existence  of  slow  travelling  waves. 

We  shall  now  obtain  slow  travelling  wave  solutions  of  (2)  which  are  in  some 
sense  near  the  stationary  heteroclinic  orbit  connecting  P  and  R  .  The 
technique  is  to  introduce  a  parameter  y  into  the  equations  which  is  propor¬ 
tional  to  E(P)  -  E(R).  In  particular,  suppose  that  H(P)  =  H(R),  where 
VH (U)  =  F (U, 0) ,  and 

F(U,y)  =  (uM (U)  +  uy  ,  vN (U) )  , 

where  F  is  as  in  case  (b) .  We  then  linearize  the  travelling  wave  equations 
about  the  above  orbit  when  y  =  0  and  apply  the  implicit  function  theorem. 
Travelling  wave  solutions  are  obtained  as  follows.  Let  W  =  H2(1R1 )2  and 

Y  =  L2(IR1)2.  Note  that  the  rest  points  of  F(U,y)  vary  with  y  .  Let  P^ 

r 

and  R^  be  the  rest  points  "P"  and  "R"  for  F(U,y)j  P_^  and  R^  depend 

,  £  £  -i 

smoothly  on  y  .  Let  A(£)  =  e’(l  +  e ')  ,  and  define 

B(£,y)  -  (1  -  A(£)) <R  -  R0)  +  MO  (P^  -  PQ)  . 

Now  let  (U,A)  be  the  solution  of  (4),  (6)  when  y  =  0,  and  let 
2 

F  :  1R  x  W  x  y  be  defined  by 

F(y,6,W)  =  (W  +  U  +  B(-,y>>"  +  F(W  +  U  +  B(.,y),y)  -  9  (W  +  U  +  B  ( • ,  y) ) 1  . 

Note  that  from  the  manner  in  which  B  was  chosen,  we  have  that  F  (W  +  B  ( • , \ )  +  (U, . ' 
e  Y  for  all  W  e  W.  Thus  F  is  a  (smooth)  mapping  of  the  indicated 
spaces.  Moreover,  by  construction  and  the  Riemann-Lebesgue  lemma,  (y,v,U) 
will  be  a  solution  of  (6)^,  (7)  if  and  only  if  U  =  W  +  B(*,y)  +  U  where 


F(y,0,W)  =  0 


Now  let  L  :  E  *  w  +  y  be  the  Frfichet  derivative  of  F  with  respect  to 
(0,W)  evaluated  at  (y,0,W)  =  0  ,  i.e., 

(14)  L(0,W)  =  W"  +  dyFOJfCDW  -  6A  , 

where  A  =  U'.  We  shall  also  need  to  consider  the  densely  defined  operator  L, 
on  y  with  domain  W  defined  by 

LQW  *  W"  +  dyFtU.OJW  . 

For  an  operator  T  ,  let  I (T)  denote  the  spectrum  of  T  . 

LEMMA  4.1.  Z (Lq)  n  {a  >_  0}  consists  of  isolated  eigenvalues  of  finite 
multiplicity. 

Proof .  We  modify  an  argument  given  by  Bardos ,  Matano,  and  Smoller,  [1],  for  a 
scalar  equation.  Let  M(5)  be  the  matrix  valued  function  d^  F(U(C),0).  Then 
M(£)  is  a  symmetric  matrix  for  all  £  ,  and  all  the  eigenvalues  of  M(±°°)  are 
strictly  negative.  Let 

Mg  ( 5)  =  X(£)  M (®)  +  (1  -  X  (0  )  M(-«)  . 

it  is  easily  verified  that  MQ(£)  is  symmetric  and  has  eigenvalues  which  are 
strictly  less  than  some  negative  constant  mQ  for  all  £  e  1R1  .  Thus  if 
rq  =  d2/d£2I  +  Mq(C),  and  if  for  U,  W  e  C^IR1)2  we  let 
ax(W,U)  =  (R0W,U)y  +  X(W,U)y  , 

then  a^  defines  a  continuous  bilinear  form  on  H^IR1)2,  and  from  the  above 

properties  of  MQ,  we  have  that 

a.  (W,W)  <  (X  +  max(-l,mn)  )!lwll2 
X  °  Hi 

Thus  a^  is  coercive  on  H^dR1)2  for  all  X  ^  o  ,  so  that  by  the  Lax-Milgram 


theorem,  RQ  +  XI  is  an  isomorphism  of  hMir1)2  onto  H'^dR1)2.  Moreover, 
if  (Rn  +  X)U  =  f  e  Y  then  U  e  W  .  Thus  Rn  +  X  maps  W  onto  Y  so  that 


-X  ^  Z(Rg).  For  such  A  we  have  that  (LQ  +  X) (rq  +  X)"1  =  I  + 


where 


=  (Mg  -  M) (Rg  +  X)  is  an  operator  on  Y  .  From  our  choice  of  Mg  we  have 


that  K.  is  compact  for  X  <_  0  >  clearly  -1  e  I(K,)  if  an<3  °nly  if  r(L.) 


The  result  follows  from  the  Riess-Schauder  theory  for  compact  operators;  see 


[1]  for  details.  D 

LEMMA  4.2.  The  kernel  of  Lq  is  one-dimensional  and  is  spanned  by  the  vector 
AU)  .  The  range  of  LQ  is  the  orthogonal  complement  of  A  iri  Y  . 

Proof.  Clearly,  A  e  ker  L  ,  so  that  if  XQ  is  the  largest  element  of  Z(Lq), 
we  have  that  XQ  >_  0.  Thus  by  Lemma  4.1,  Xq  is  an  element  of  the  point  spec¬ 
trum,  so  that  if  Q (W,W)  =  (LpW.W)^,  then 

XQ  =  sup{Q(W,W)  :  llwll  =  1,  W  £  W}  , 

and  this  supremum  is  attained  at  a  vector  WQ  e  W.  Suppose  that  XQ  >  0  .  Then 
we  must  have  that  (Wq,a)y  =  0.  However,  if  has  components  a  and  x  , 

then 

(15)  XQ  =  /  (-a2  -  t2  +  fua2  +  2fvat  +  gvT2)d^, 

—  <30 

where  the  partials  f  ,  f  ,  and  g  are  evaluated  at  U (2) .  Now  f  (u,v)  <  0 

U  V  V  '  v  — 

for  all  u,v  >_  0 .  If  o  and  t  ever  had  the  same  sign  and  were  nonzero  for  some  £  , 
we  would  get  a  strictly  larger  value  in  the  integral  in  (15)  by  replacing  a 
with  | a  I  and  t  with  -|t|.  Thus  a  and  x  must  always  have  different  signs. 
Moreover,  neither  a  nor  t  ever  change  sign,  since  if  this  was  the  case  at 
x  =  ,  ( I o [ , —  |  t ( )  would  again  be  an  eigenvector  of  LQ.  However,  by  the  above, 

a(£Q)  =  t(£0)  =  0  ;  in  addition,  o'  (CQ)  =  T,(£Q)  =  0  ,  since  |o|  and  — [ x | 
must  both  be  locally  H2  functions.  Since  L  (|o|,-|x|)  =0  is  a  homogeneous 
second  order  system  of  o.d.e.'s,  these  conditions  imply  that  o  and  x  vanish 
identically.  Thus  a  and  x  must  be  of  constant  and  different  sign.  However, 
by  Lemma  3.4,  a  similar  remark  holds  for  the  components  of  A  .  Thus  (A,WQ)y 
must  be  nonzero,  so  that  XQ  =0.  A  similar  argument  shows  that  zero  is  a  simple 
eigenvalue  of  LQ. 


22' 


The  last  assertion  of  the  lemma  follows  from  the  fact  that  range  (LQ)  = 
range  (L^k  , ')  is  closed  in  Y  ,  and  the  usual  relation.  A1  =  (ker  L*)1  = 
cl (range (LQ) ) .  □ 

Now  let  W  =  A1,  and  define  a  mapping  fjIR^  x  ■>  y  by 
f(y,Tr,0,W1)  =  F(y,0,ttA  +  Wj). 

THEOREM  4.3,  There  exists  a  smooth  function  h:lR2  -»•  IR  x  ^  where  h(0)  =  0, 
and 

h(y,ir)  =  (6  (y,tt)  ,Wj  (y,h)  )  , 

such  that  the  set  of  solutions  of  f (y,it, 0,Wj)  =  0  near  the  origin  consists 
precisely  of  vectors  of  the  form  (y,h-,0  (y,u)  ,Wj  (y,tt)  ) . 

Remark.  The  additional  parameter  u  is  present  due  to  the  inherent  degeneracy 
in  the  problem  arising  from  the  fact  that  translations  of  solutions  are  solu¬ 
tions.  we  can  therefore  "normalize"  the  space  by  setting  u  =  0  ,  and  by  the 
uniqueness  assertion  of  Theorem  4.4  we  are  assured  of  obtaining  all  solutions 
near  the  origin  modulo  translations. 

Proof.  We  have  that  the  derivative  of  f  with  respect  to  (OjWj)  is  precisely 
the  restriction  L  of  L  (defined  in  (14))  to  IR  *  Il/j  .  The  theorem  will 
follow  from  the  implicit  function  theorem  if  it  can  be  shown  that  Lj  :  JR  *  ■*  Y 

is  an  isomorphism.  If  Li(0,Wj)  =  0,  we  find  by  multiplyinq  this 
equation  by  A  and  integrating  that  -9(A,A)y  =  0  ,  so  that  9=0,  and  that 
W  <■  ker  Lq  n  Wj  =  {0}.  Hence  L ^  is  1:1.  Now  let  G  e  Y  ,  and  put  9  = 

- (G,A) y/(A,A)y  .  Then  1^(0, W  )  =  G  has  a  solution  if  and  only  if 
(16)  LQW  =  G  +  9A 

has  a  solution  in  C i/  .  From  our  choice  of  0  and  Lemma  4.2,  this  is  obviously 
the  case;  thus  L^  is  onto  .  O 

COROLLARY  4.4.  0  (0,0)  <  0  . 

Y 


23- 


-V*  "i- 


»V*-, 


Proof 


Differentiate  F<y,6 (y,0) ,Wj (y.O) )  =  0  with  respect  to  y  to  obtain 

r"  +  Mr  +  f  (u,o)  -  0  (o ,oj a  =  o  , 
y  y 

where  r  =  W1v(0,0)  +  B  (• ,0) .  Note  that  F  (U,0)  =  (u,0),  where  u  is  the 
IT  Y  Y 

first  component  of  (J  .  Multiply  the  above  equation  by  A  and  integrate  to 
obtain 

/°°  A  •  (u,0)d£  -  0  (0,0)  /°°jAj2d£  =  0  . 

_oo  — oo 

Thus ,  we  have  that 


n^foo)  —  ii ^  oo)  _ 

0  (0,0)  =  - — — -  <  0  .  c 

Y  2||  All  2 

We  remark  that  all  of  the  above  results  apply  also  to  the  mapping  F  3  = 

F  +  y2G(VI)  where  Gs  W  -*■  Y  is  a  smooth  mapping.  Hence  travelling  wave 

solutions  persist  under  small,  smooth  perturbations  of  the  nonlinear  term.  In 


particular,  we  see  that  suv.h  behavior  is  not  peculiar  to  Hamiltonian  systems. 


5.  Existence  of  stable  dines. 

In  this  section,  we  introduce  a  spatially  inhomogeneous  perturbation  into 
the  equations.  It  is  no  longer  true  that  translations  of  solutions  are  solu¬ 
tions;  indeed,  for  a  reasonable  class  of  perturbations,  the  spectrum  of  the 
linearized  perturbed  equations  moves  into  the  stable  half  plane.  Thus  even  a 
slight  inhomogeneity  in  the  relative  competitive  strengths  of  the  two  species 
can  cause  the  densities  to  evolve  to  a  stable  large  amplitude  pattern.  This 
fact  has  been  observed  by  Peletier,  [16],  and  by  Peletier  and  Fife,  [17],  for  a 
scalar  equation  arising  in  population  genetics.  Mimura  and  Nishiura,  [14], 
have  obtained  small  amplitude  spatial  patterns  for  a  system  arising  in  develop¬ 
mental  biology  by  introducing  spatial  inhomogenities  via  perturbed  bifurcation 


theory . 


We  will  assume  that  H(P)  =  H (R) ,  where  VH(U)  =  F(U,0)  and  that  f  , 

(case  b) ,  has  the  form 

(17)  F(U,x,e)  =  (uM(U)  +  (ePj  +  e^jJUjVN  (U)  +  (ep2  +  e2p2)v)  ; 
here  p^,<]k  e  L2(1R1),  i  =  1,2.  We  shall  apply  an  argument  based  on  the 
Liapunov-Schmidt  procedure  to  obtain  solutions  of  (1)  which  also  satisfy  (5). 

We  shall  need  the  following  theorem. 

THEOREM  5.1.  Suppose  that  G  is  a  smooth  mapping  of  a  neighborhood  of  the 
origin  in  1R  x  W  into  Y  ,  where  CU  and  Y  are  Banach  spaces.  Suppose  also 
that  G(0,0)  =  0  ,  and  for  (s,W)  e  ]R  *  W  that 

(i)  G__(0.0)H„  =  -G  (0,0)  for  some  H„  e  W  , 

w  o  e  - o 

(ii)  ker  G..(0,0)  «  span  {a},  for  some  A  ^  0  in  W  , 

w  - 

(iii)  the  range  of  G^(0,0)  is  closed  in  Y  and  has  codimension  one, 

(iv)  there  exists  a  one  parameter  family  W^(s)  e  W  for  small  s  with 
WjtO)  =  W^s  (0)  *  0  such  that  G(0,sa  +  W^s))  =  0  , 

(v)  6eW(0,0)A  +  6ww<0'0> <H0,A>  ^  range  (GW(0'0))* 

Then  the  set  of  solutions  of  G(e,W)  =  0  near  the  origin  in  IR  x  W  consists 

precisely  of  two  transversally  intersecting  curves.  The  nontrivial  branch  has 
the  tangent  (1,H),  where  H  =  Hq  +  YA  for  some  y  e  R1  so  that  it  can  be 
parameterized  by  e  .  We  therefore  obtain  a  curve  W(e)  e  W  such  that  W(0)  =  0 

w  (0)  =  -H,  and  G(e,W(e))  =  0  . 

e 

The  proof  of  a  similar  theorem  can  be  found  in  Nirenberg,  [15,  Theorem  3.2.2] 
we  shall  therefore  only  sketch  the  details.  Conditions  (i)  -  (iii)  and  the 
Liapunov-Schmidt  procedure  can  be  used  to  show  that  solving  G(e,w)  =  0  near 
the  origin  in  ®  *  111  is  equivalent  to  solving  G(a,b)  =  0  near  the  origin 
in  1R2 ,  where  G(a,b)  *  y*G(a,aHQ  +  bA) ,  and  y*  is  a  bounded  linear  func¬ 
tional  on  Y  which  annihilates  the  range  of  0,0).  Since  VG(0,0)  *  (0,0), 
the  Morse  lemma  can  be  used  to  obtain  the  desired  result  provided  that 


d^G (0 ,0)  = 


gI  g2 
g  2  g3 


is  nondegenerate  and  indefinite,  i.e.,  that  g  ^g  3  -  g2  <  0. 


Hypothesis  (v)  implies  that  g2  ^  0  whereas  (iv)  implies  that 
(18)  g  =  G..  (0,0)  =  y* (G  (0,0) (A, A) ) 


-  (y*G(0,sA  +  W. (s) ) 

,  ■>  1  s=0 


Finally,  the  tangent  of  the  trivial  branch  is  (0,A);  in  addition,  the  tangent 
of  the  nontrivial  branch  at  the  origin  is  transverse  to  that  of  the  trivial 
branch  and  lies  in  the  plane  spanned  by  (1,HQ)  and  (0,A);  (see  Nirenberg) . 
Hence  the  nontrivial  branch  can  be  parameterized  by  e  in  the  manner  described 
above. 

How  define  a  (smooth)  mapping  (3:1  x  111  +  V  by 
G(e,w)  =  (U  +  w)"  +  F (U  +  w,  x,  e)  , 
where  F  is  as  in  (17). 

THEOREM  5.2.  Suppose  that  Pj  and  are  chosen  such  that 

(19)  /  (a(x)u(x) p1 (x)  +  b(x)v(x) p2(x) )dx  =  0. 

—  00 

Then  for  an  open  set  of  (p^p^)  e  Y  lying  in  the  hyperplane  defined  by  (19) 

there  exists  a  solution  W(e)  of  G(e,W)  =  0  ,  for  sufficiently  small  e  . 

Proof.  We  show  that  the  hypotheses  of  Theorem  5.1  are  satisfied.  From  Lemma 

4.2,  we  see  that  (19)  implies  that  G£(0,0)  =  (uPj,vp2)  is  in  the  range  of 

G  (0,0),  i.e.,  A1;  (ii)  and  (iii)  follow  from  the  same  lemma.  Condition  (iv) 

w 

is  verified  by  noting  that 

-2  x+s 

U(x+s)  -  U(x)  =  s  A  (x)  +  s2{s  /  jZ-  F(U(n>  )dnde)  ; 

X  X 

the  expression  in  brackets  tends  to  -F(U(x))  uniformly  on  compact  sets  as 
s  ■+•  0  .  Moreover,  the  derivatives  of  U(x)  are  uniformly  bounded,  and  tend 
to  zero  exponentially  as  )x|  -+  00 .  It  easily  follows  that  the  expression  in 
brackets  is  uniformly  bounded  in  W  as  s  ->•  0  . 


-26. 


We  finally  show  that  (v)  holds  for  most  (p  ,p  )  that  satisfy  (19). 

1  2 

Note  that  if  Ln  =  G^(0,0),  then  L^  is  a  bounded  self-adjoint  linear  operator 

depends  linearly  and  continuously  on  (p^p^). 


on  A1  so  that  HQ  =  Lq1 


up 


vp 


2  > 


Condition  (v)  is  equivalent  to 

(20)  0/  /“< a2p  +  b2P2)dx  +  /°°Atd2F(U(x) ,0)  (A,H0)dx. 

—00  —00 

Let  V  =  d2F  (U (x) ,0) (A,A) ;  from  (18),  we  see  that  V  £  A1  ;  hence,  using  the  fact 
that  L  1  is  a  self-adjoint  operator  on  A1,  we  have  that 


(v'Vy  = 


UP] 

vp. 


Moreover,  we  have  that  if  U(x,s)  =  LJ(x+s),  then 


ds 


2  (Uxx  +  F (U  ,0) 1 


=  0  , 


s=0 


so  that  if  G  =  U  (x,0)  =  -F(U(x)),  we  have  that 
ss 

LqG  +  d2F (U, 0) (A, A)  =  0  ; 

that  is,  G  =  -L^v.  Thus  L  *V  =  P(U),  and  (20)  is  equivalent  to 


(21) 


0  ^  /  [(a2  +  uf(U))p  +  (b2  +  vg(U))p2)dx  . 


This  fails  to  hold  for  all  (p  ,p2)  which  satisfy  (19)  if  and  only  if 
(a2  +  uf,  b2  +  vg)  =  s(ua,vb) 

for  some  s  £  IR .  However,  in  this  case  at  least  one  of  a2  -  sua,  b2  -  svb 
must  be  of  constant  sign,  whereas  -uf(U)  and  -vg(U)  both  change  sign.  □ 
Remark.  Since  the  second  order  perturbations  <p1  and  <t>2  are  free  parameters, 
we  see  that  Theorem  5.2  actually  holds  for  a  generic  class  of  L2  perturbations. 
COROLLARY  5.3.  Let  X(t)  be  the  largest  eigenvalue  of  the  problem 
(22)  Z"  +  dyF(U  +  W(e) ,x, e) Z  =  XZ, 


and  let  Z  =  1(e)  be  a  corresponding  eigenvector.  Then  for  all  (p1#P2) 
satisfying  (19)  and  (21),  we  have  that  A£(0)  ¥  0  . 


Proof. 


We  have  that  W  (0)  =  H„  +  yA  for  some  y  e  3R.  Let  Z  =  Z  (0) 

e  0  e  e 

and  A  =  A  (0) .  Differentiate  (22)  with  respect  to  e  ,  set  e  =  0  ,  multiply 

e  e 

by  A  ,  and  integrate  to  obtain 

2  oo  f 

A£  =  II  Ally  [  f° I (a2P1  +  b2p2)dx  +  /  A  d2F(U,0)(W  (0),A)dx); 

— OO  -00 

we  have  used  the  fact  that  A (0)  =  0  and  that  Z (0)  =  A  .  From  (18)  and  the 
fact  that  W^(0)  =  HQ  +  yA  ,  it  can  be  seen  that  the  expression  in  brackets 
equals  the  right  hand  side  of  (21).  □ 

Thus  for  e  of  the  appropriate  sign,  the  spectrum  of  the  linearization  of 
the  perturbed  equations  will  lie  in  the  stable  half  plane.  Since  the  linearized 
operator  is  clearly  sectorial,  standard  abstract  results  can  be  invoked  to 
obtain  the  full  nonlinear  stability  of  the  solutions  of  the  perturbed  equations 
with  respect  to  appropriate  perturbations  in  the  initial  data  of  (2)?  see  for 
example,  Henry,  [8],  Indeed,  such  solutions  must  be  asymptotically  exponentially 
stable.  We  also  note  that  such  stable  behavior  persists  under  small  c1  per¬ 
turbations  in  the  nonlinear  term;  this  is  an  immediate  consequence  of  the  implicit 
function  theorem.  Even  though  such  perturbations  may  have  to  be  extremely  small, 
we  nevertheless  obtain  the  structural  stability  of  stationary  solutions  connecting 
P  and  R  .  This  fact  was  not  clear  from  our  discussion  of  autonomous  systems 
in  section  3. 

6.  Concluding  Remarks . 

Our  results  indicate  that  diffusion  coupled  with  a  nonlinear  competitive 
interaction  can  lead  to  the  formation  of  large  amplitude  patterns  which  are 
locally  consistent  with  the  principle  of  competitive  exclusion  and  yet  on  a 
global  scale  allow  both  species  to  persist.  (Although  such  results  were  proved 
under  the  hypothesis  of  gradient  or  near-gradient  structure,  the  theorems  of 
section  2,  (proved  in  the  absence  of  such  hypotheses),  suggest  that  such  behavior 
is  exhibited  by  solutions  of  quite  robust  classes  of  equations).  Furthermore, 
stability  of  stationary  patterns  was  achieved  through  the  introduction  of  spatial 
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inhomogeneity  in  the  relative  competitive  strengths  of  the  two  species.  It  is 
likely  that  this  is  also  a  necessary  condition  for  stability  of  such  patterns. 

In  particular,  the  technique  of  Bardos  and  Smoller,  [2],  can  probably  be 
applied  to  obtain  the  instability  of  the  homoclinic  orbits  and  the  solutions 
of  the  Neumann  problem  obtained  for  autonomous  systems  in  section  3.  It  would 
be  interesting  to  see  whether  stable  solutions  of  the  singularly  perturbed 
Neumann  problem,  (9),  can  be  obtained  through  the  introduction  of  a  second, 
spatially  inhomogeneous  perturbation. 

Another  interesting  question  is  the  introduction  of  time-dependent  perturba 
tions.  Forced  oscillations  in  competition  equations  without  diffusion  were 
studied  numerically  by  Koch,  [101.  It  therefore  seems  natural  to  introduce 
time-dependent  periodic  perturbations  into  systems  that  admit  stable  large  ampli 
tude  stationary  solutions.  Slowly  varying  perturbations  are  likely  to  give  rise 
to  stable  large  amplitude  periodic  solutions,  whereas  rapidly  varying  perturba¬ 
tions  may  give  rise  to  more  complicated  recurrent  behavior. 
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